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On the night of May 23/24, 1975, a sequence of rocktt and balloon experiments was launched 
from Chilca Base, Peru (12.5°S, 76.8°W, magnetic dip - -0.7°) as part of Operation Antarqui. Detailed 
analysis and comparisons of the data have yielded a significant and unexpected result , namely the first 
direct measurement of lower mesospheric response to a galactic x-ray source. This result could only have 
been determined at the equator, where cosmic ray background effects arc minimal. 

Our objective was to seek out the equatorial energetic electron belt, sporadically reported to con- 
tain fluxes near auroral levels (e.g., Heikkila. 1971 ; Goldberg, 1974), measure the bremsstrahlung radiation 
produced by this particle belt, and determine the influence of this radiation on the middle atmosphere. It 
was assumed that if the quasi-trapped particle belt were present, it would be subject to the nighttime downward 
drift forces which lower the equatorial F region ionosphere. The particle belt wouid therefore drop into a 
more dense region of the atmosphere where significant enhancement of bremsstrahlung radiation production 
would occur. Two high altitude rocket payloads (Nike Tomahawk 18.170 and 18.171) were launched to 
probe the thermosphere during and following the anticipated downward drift period. Each carried an on-axis 
x-ray scintillation detector with a 5 mil Be window for energetic pat tide rejection, and with pulse height 
discrimination in the energy ranges 5- 1 0 (X 1 ). 10-20 (X2), 2040 (X3L and >40 keV (X4). Geiger Mueller 
energetic electron detectors with a multiplicity of window thicknesses and look directions were carried to 
monitor energetic electron fluxes above thresholds between 15 and 40 keV as a function of pitch angle. 
Magnetometers and lunar sensors were used to determine payload aspect. 

Each Nike Tomahawk was accompanied by racket probes to sample conductivity (o) and ozone 
(Oj) below 80 km, to search for modifications in these parameters caused by the thcrmospherically produced 
x-rays. Meteorological datasondes (MET) were also used to evaluate the meteorological characteristics of the 
middle atmosphere and provide control data during the period of study. Finally, monitoring of cosmic rays 
and energetic x-rays near 40 km was made on a continuous basis during the entire night with balloon-borne 
instrumentation. The precise sequencing of each measurement is depicted in Figure I. 

Figure 2 compares the count rates for each of the x-ray channels during the flights of 18.170 and 
18.171. The vertical line displayed at each 10 second interval represents the extremes of coun* rate for each 
second determined in the interval t ± 5 seconds. Midpoints of each range are connected to permit tracking 
of the mean flux through the fligh* Both time and altitude arc displayed for reference. The sharp rise in 
count rate near 50 seconds represents the sudden exposure of the x-ray detector to ambient radiation upon 
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nose cone deployment. Finally, the count rate enhancement on X4 lot each flight above 3U0 km is attributed 
to energetic electrons (>150 keV) penetrating the Ik' windows This flux is caused by the proximity of the 
launch site to the outer fringe of the South Atlantic anomaly. 

The low energy channels show a significant enhancement in count late between the first and second 
measurement of x-ray flux. It reaches a maximum ratio of 4.5 on XI . and is not present on X4. suggesting 
the added contribution to be caused by a relatively “soft” source. The enhancement has now been analyzed 
in sufficient detail to permit full accountability by x-ray galactic sources, with no need for a bremsstrahlung 
contribution. The absence of an equatorial energetic electron belt during this period, at the minimum of the 
sunspot cycle, has been further verified by near zero count rates obtained with all electron detectors aboaid 
both flights, although each (light remained above 300 km for more than 150 seconds. Woodman (1975 ) has 
also indicated that the major portion of downward drift occurred prior to the 18.170 launch, using Jicamarca 
Observatory drift data. s would imply that the particle belt (if present) would have been near its lowest 
altitude during both me- .icments, thereby improving the chance of payload penetration. Finally, the 
spectral distribution and intensity of the x-ray flux observed aboard I8.I7C has been analyzed and found 
to agiee favorably with the diffuse galactic x-ray backgr jund. 

The enhancement observed between flights can he attributed to galactic stellar x-ray sources, 
dominated by the presence of Sco X-l , which was at a suitable location for detection and is 10 times brighter 
than any other x-ray star in the 5-10 keV range. Figure 3 is a polar plot using elevation and azimuth co- 
ordinates relative to the 18.171 payload position. We have depicted the look direction of the payload (and 
x-.ay detector), and the positions of the moon , Sco X-l . and galactic center in this reference system. 

Although the trajectory azimuth direction is nearly d e west, an irregularity during nose cone deployment 
caused the payload to tip in the backward direction, and stabilize with a ±5° cone centered near 125° azimuth 
and 60° elevation. As will be seen, this was highly fortunate, as it permitted isolation of Sco X-l as the 
primary source responsible for the observed x-ray enhancement. 

A power spectrum analysis, using fast Fourier transform techniques, was applied to the data of each 
x-ray channel on both flights. Because of the backward tip of the scintillator, the angle between detector 
and Sco X-l was 46 ± 5°, which is sufficiently large to detect modulation effects in the lower energy channels 
ductoconing. Figure 4 shows the spectral plot obtained from the XI data on 18.171. A singular peak at 
k index (sampling interval/pcriod) = 100 at 9 db above the RMS noise is apparent in the plot, and relates 
precisely to the 2.27 second coning period of the payload Here, a unit shift in k is equivalent to about a 1% 
change in period. Other channel on this flight exhibited the same result, but with less definition No 
periodicities could be discerned in the data of 18.170. The amplitude spike denotes the net contribution 
from one or more point sources ( A.C. component). Further analy sis of phasing Jata lias established the 
point source to be within 10° azimuth of Sco X-l . Finally, in Figure 5 we show the energy spectrum for 
the A.C. component of x-ray flux measured on 18.171 . and compare it with the average spectrum for Sco X-l 
(Laros and Singer, 1976). The observed spectrum is found to be slightly harder than that of Sco X-l and 
somewhat more intense. However, Holt et. al. (1976) have recently reported wide variability in the 3-6 keV 
intensity flux emitted front Sco X-l for a four month pcru>d just prior to this measurement, using data from 
the Anel-5 All Sky Monitor. This range is indicated on the figure to demonstrate that our extrapolated 
spectrum lies within the reported unge. Finally, the absence of periodic data from 18.170 i- consistent 
with Sco X-l being low in the eastern sky and out of detector range at that time. 

We now consider the stratospheric measurements which demonstrate an atmospheric response to 
the measured Sco X-l radiation. Figure 6 shows the two nighttime conductivity profiles (2000 and 0153 
LMT) obtained in parallel w.th die high altitude measurements. In addition, a daytime value on May 28 
(1423 LMT) is illustrated to demonstrate measurement consistency Below 55 km, all three profiles are 
nearly identical, and reflect the nearly constant influence of equatorial cosmic rays as an ionization source 



in this region. The early nighttime profile continues to track the daytime value up to 65 km, above which 
solar Lyman a causes the daytime conductivity to grow with altitude at a far mure rapid rate. Our primary 
concern is the enhancement of the late nighttime profile between 55 and 75 km, causing a factor of two 
increase near 65 km. This increase is attributed to the Sco X-l contribution obser ed aboard 18.171. fhe 
height of energy deposition is consistent with the calculations of Berger et. al. (1974), which predict that 
5-10 keV x-rays will be subject to maximum atmospheric absorption between 60-75 km. Finally, we note 
that the measurements of meteorological parameters and ozone during the period of conductivity measure- 
ment do not indicate dynamical behavior which could account for the observed conductivity enhancement. 

We next turn to aeronomic considerations to establish consistency with the preceding interpretation. 
Th e ii ir.ie dependent continuity equation near 60 km may be written as 

(dN/dt)*= Q -a t N- (1) 

where N is total plasma density, a 4 is the ion-ion recombination coefficient, and Q is the ionization source 
function, composed of an x-ray (Q X r) and 3 cosmic ray (Q CR ) component. The steady state solution of 
this equation is 

N 2 * Q/o* (2) 

Assuming that the first measurement (t 0 ) was made near equilibrium for Qcr and the second (t 2 ) near 
equilibrium for Q CR + Oxr. we obtain 


(N 2 /N 0 )2 = 1+(Ox R /Ocr) (3) 

where N 2 is the steady state value at t 2 and N 0 , the value at t 0 . Since N 2 /N 0 = 2, Q X r/Qcr 3. The 
calculated Q XR (65 km) = .006 ion pairs cm' 3 s _ > using the date of 18.171 . At first glance, the extrapolated 
Q cr from the cosmic ray ionization chamber aboard the high altitude balloon appears to be of equivalent 
magnitude. However, Williamson (1976) has shown that for lower latitudes, the cosmic ray ion production 
rate exhibits a more rapid decrease with increasing altitude than has been assumed in the past. This empirical 
result is based on published data from numerous balloon flights, and implies that a value for Q CR (65 km) 

~ .001 would be more appropriate. 

The desirability for a low Q CR near 65 km is even more acute when the full solution of (1) is taken 
into account. Under the assumption that Q and a 4 remain constant between flights, we have 

N = C(A - exp (-Bt))/(A + exp (- Bt)) (4) 

where 

CMQ/aj), B = 2(Qa i ) ,/2 . 

A * (C -N 0 )/4C + N 0 ). 

•i d N 0 is the initial value at t = 0. We have conducted a parametric study of this solution to evaluate the 
conditions present during the late evening measurement (0153 LMT, t = 6 hrs) assuming that t = 0 represents 
the early measurement (2000 LMT). For these calculations, N 0 is 150 cm -3 as estimated from the blunt 
probe data, and a is of order 5 X 10‘ 8 cm 3 s' 1 . The results are shown in Table 1 , where we find that N, 

6 (hrs) represents from 63 to 68% of the equilibrium value. Since the requirement that N () be doubled can 
be met during this interval, this implies that at steady state, N 2 /N 0 be larger than 2, and 0 ( - R be sufficiently 
small to meet the criterion. The values in the table are reasonable within the experimental and theoretical 
guidelines of this work. Decreasing a or increasing Q will also permit the doubling effect to occur, but by 
exceeding limits considered reasonable for the appropriate parameters. 
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Table I 
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(cm 3 s'* ) 

0 

(cm' 3 s' 3 ) 

No 

(cm' 3 ) 

N, (6 hrs) 
(cm -3 ) 
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(cm' 3 ) 

% 

(100X n,/n 2 ) 

5 X I0- 8 

7.5 X I0* 3 




68 

4 

7.5 

■ ... 

mam, 


63 

5 

8.0 

150 

271 

400 

68 

5 

10.0 

150 

305 

447 

63 

4 

9.0 

150 

298 

474 

63 


Finally, other experimental results have claimed observed modulation of radio waves caused by 
Sco X-l nighttime enhancement of upper D region ionization (e.g., Edwards et. al., 1969; Ananthakrishnan 
and Ramanathan, 1969; Chilton and Crary, 1971) using limited statistical correlations. These results in 
turn have been discussed and disputed on theoretical grounds by Poppoff and Whitten (1969), Francey 
(1970), and Poppoff et. al. ( 1975). Although our results pertain to modifications of ion density at altitudes 
below which electrons can usually survive at night, it would appear that these new results show sufficient 
cause to reconsider the earlier theoretical arguments regarding ionospheric response to x-ray galactic sources. 
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Figure 1. Sequencing of Measurements 
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